Effects of interactions in transport through Aharonov-Bohm-Casher interferometers 
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We study the conductance through a ring described by the Hubbard model (such as an array of 
quantum dots), threaded by a magnetic flux and subject to Rashba spin-orbit coupling (SOC). We 
develop a formalism that is able to describe the interference effects as well as the Kondo effect when 
the number of electrons in the ring is odd. In the Kondo regime, the SOC reduces the conductance 
from the unitary limit, and in combination with the magnetic flux, the device acts as a spin polarizer. 
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Advances in semiconductor technology have provided 
useful tools to test fundamental concepts of quantum 
physics, such as the superposition principle and the exis- 
tence of topological phases [l[ . Beautiful demonstrations 
of these are studies of the Aharonov-Bohm (AB) effect 
in mesoscopic rings, particularly with embedded quan- 
tum dots (QDs) 0, [j]. The effect of interactions in these 
systems is still a matter of debate [Hj. Despite the enor- 
mous effort to describe transport through interacting re- 
gions 0], at present we do not have a unified procedure 
to extend the results of the single particle case to many- 
body cases. A serious shortcoming of the calculation of 
the conductance G through an interacting ring is that 
even knowing the exact eigenstates of the ring, there is 
no simple procedure to calculate G. When the coupling 
V of the ring to the conducting leads is small, Jagla and 
Balseiro (JB) used a perturbative expression in V for 
G that is exact for any V in the non-interacting limit 
0. Similar equations were used recently, assuming that 
a Zeeman term destroys the Kondo effect in the system 
[1, iSi] . Another expression in order was proposed last 
year Unfortunately these expressions are not valid 
in the Kondo regime, in which the number of electrons 
in the ring is odd, because the resulting Kondo physics 
cannot be described by perturbation theory in V. The 
ideal conductance in the Kondo regime was recovered by 
mapping the model into an impurity Anderson model, 
but in this formulation interference effects were lost 0]. 
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Recently, the Aharonov-Casher (AC) effect 
charge-spin dual of the AB effect, has been demonstrated 
exp erimentally in semiconductor mesoscopic rings |12l . 
|l3|. The AC phases are originated due to the Rashba 
spin-orbit coupling (SOC) in the ring, resulting from 
electronic motion in the presence of an electric field nor- 
mal to the plane of the ring. The interference between 
electrons of given spin travelling clockwise and anticlock- 
wise produces a strong modulation of the electronic cur- 
rent through the device. Recent theoretical research 
has successfully explained the modulation of the conduc- 
tance in terms of non-interacting electrons. However, 
the single-electron picture, turns out to be inadequate to 
describe electronic transport in the strongly-interacting 



case, particularly in the Kondo regime, as we will show. 

In this Letter, we describe a systematic procedure to 
calculate equilibrium conductance G through a ring of an 
interacting system weakly coupled to conducting leads, 
that takes into account both the effects of interference 
and correlations in presence of a magnetic flux and SOC. 
Using a non-abelian gauge transformation (NAGT), we 
show that for on-site interactions, the SOC can be ab- 
sorbed in opposite AC phases for spin up and down in 
an adequately chosen quantization axis. For a Hubbard 
model (that describes a ring of an even number of QDs), 
in absence of SOC G vanishes when the magnetic flux 
amounts to half a flux quantum. For other fluxes in the 
Kondo regime, G reaches the unitary limit (ideal con- 
ductance [J]). When the SOC is turned on, the ideal 
conductance is destroyed and G shows a strong spin de- 
pendence in this regime. 

Our first task is to derive the appropriate extension 
to the Hubbard model to include the SOC in an ade- 
quate representation that simplifies our subsequent cal- 
culations. To illustrate the procedure, it is easier to be- 
gin with non-interacting electrons in the continuum. The 
correct Hamiltonian for this case was derived by Meijer et 
al. [3. The SOC is Hsoc = a'a-IS x (p-eA), where a' 
is the Rashba constant and E is the electric field, which in 
our case is in the z direction, perpendicular to the plane 
of the ring. Including SOC and the orbital effects of the 
magnetic field, but neglecting the Zeeman term (usually 
several orders of magnitude smaller than the Kondo en- 
ergy scale in QDs Q), the Hamiltonian can be written in 
the form Q (a) 

n 2 
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where = h/ (2m*r ), m* is the effective electron mass, 
r is the radius of the ring, ujso = oL/hr, a = ha'E^, 
(f) = Birr^ is the magnetic fiux, i/io = h/e is the flux 
quantum and ct,. (<^) = ctx cos ip + ay sin ip is the Pauli 
matrix in the radial direction, and ip is the azimuthal 
angle (see Fig. 1). Although the Schrodinger equation 
Hmxi'P) — E^xiv) (where x is a spinor) has been solved 
[l4| , we are interested in a simplification of this equation 



that can be extended to the interacting case. This can 
be achieved by a NAGT xi'P) = U{f)x'{'p), where the 
operator U (</?) satisfies the differential equation 
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It can be easily checked that in the transformed Hamil- 
tonian, H^^ = (/''HyuU = —hfld'^/ip'^ the magnetic flux 
and the SOC disappeared, and enter now in the bound- 
ary condition, since x(27r) — x(0) iniphcs x'(27r) = 
f/+(27r)x'(0). The solution of Eq.® with U{0) = 1 is 



U{ip) = exp 
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where ng ~ (— sinS, 0, cos 0), 9 = arctan (ojso/fi), and 

if' = ip^l + (Wso/fl)2. 

To construct the tight binding version of iJ^i' 
assume that we have N sites, lattice parameter a (with 
Na = 27rr) and site at angle (/? = 0. For simplicity we 
consider only hopping between nearest neighbors (NN). 
Then, we can take a constant hopping t between all NN, 
except between sites iV — 1 and 0, in which the bound- 
ary condition should be included. The matrix U (27r) is 
easily diagonalized in the quantization axis ng and its 
eigenvalues are exp[i($^B -I- a<^Ac)\, where a ~ I (-1) 
for spin up (down) in this direction, '^ab = 27r0/0Oi a-nd 
'^Ac = + {'^so/^YV^'^ - !}• Therefore, destroying 
a particle with spin a at site — 1 and creating it at 
site should be accompanied by the corresponding ex- 
ponential factors. On-site interactions arc not affected by 
the NAGT. Changing the phases of the second quantiza- 
tion operators to write the Hamiltonian in rotationally 
invariant form, the transformed Hubbard model in the 
ring becomes 
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From the curvature of the dispersion relation at small 
wave vector t = fi^ / {2m* a?) , and then lOso/^ ~ 
aN/{27rta). Thus, the AC phase can be written as 
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Therefore, for large a or N, the properties of the system 
are periodic with a as observed experimentally 12, Il3l | . 



The fact that the SOC can be gauged away in one di- 
mension has been noted previously [l6|, but the explicit 
form of the transformation has not been derived. This 



transformation has important consequences. In the ther- 
modynamic limit the boundary conditions are irrelevant 
and therefore the thermodynamic properties of the sys- 
tem should be identical to those of the Hubbard model 
without SOC. This is not obvious in alternative treat- 
ments In particular it seems that the opening of a 
spin gap in the system requires long-range interactions. 
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Figure 1: Scheme of the system. 

To study the conductance, we must consider the 
Hamiltonian of the complete system H — Hi + H' 4- 
Hv, where with the appropriate quantization axis [l8[ 
and choice of phases Hi = tc{J27=^,a + 
X^i^i cr '^1-1-1 o-CitT + H.c), describes the non-interacting 
leads, and Hy = V{J2a 4)adN/2,a 
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the coupling between the ring and leads. For simplicity 
we will focus here on the particular case = 4, illus- 
trated in Fig. [TJ We assume that the leads are described 
by a constant density of states po = l/W, and we take 
for the band width of the leads W = 60i {W is usually 
much larger than t in QD arrays). The Fermi level is set 
at ep ~ 0. To control the charge in the ring wc add to 
H^ a term ~Vg J2ia ^L'^*'^ ^^^^^ represents the effect of a 
gate voltage. Our approximations to calculate G amount 
to a truncation of the Hilbert space of H'. and a slave 
boson mean-fleld approximation for the resulting gener- 
alized Anderson model (GAM). H^ can be diagonalized 
exactly (numerically for not too large A^). We retain only 
two neighboring charge configurations with n and n — 1 
particles and we have chosen n = 4. Furthermore, we 
retain only the lowest lying singlet state for 4 particles 
(iV'o) with energy Eq) and all doublets for 3 p articles. 
This procedure is valid for small enough V hj. Calcu- 
lating the matrix elements of Hy in the truncated Hilbert 
space leads to a GAM 

Hgam ^ Hi+J2 ElJi,^,){iP^J + i?o'|Vo')(^ol + 

+V E KJV'o')('/'f.|c,. + H.c.), (6) 

J,cr,?7— 0,1 

where IV'fcr) and Ej^ denote the j-th eigenvector and 
eigenvalue of H^ in the configuration with 3 particles with 
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spin projection a, in ascending order of energy and 



Hqam can be expressed exactly in terms of a slave-boson 
representation similar to that proposed by Coleman [2^ : 

\^%){^U - 4/.- l*o)(*^l - b^b, \^%){^t\ - fib 
and |^'o}(^jCTl ~^ b'^fjcr, where the operators b'^ and 
create a boson and a fermion respectively and are subject 
to the constraint J2ja fjafjo- + b''b = 1, which is incorpo- 
rated in the Hamiltonian with a Lagrange multiplier A. 
We perform a saddle-point approximation in the bosonic 
degrees of freedom, which reproduces the Kondo physics 
at low temperatures and becomes exact in the limit of 
infinite degeneracy of the magnetic configuration, due to 
vanishing fluctuations around the mean-field value of the 
bosonic field The problem becomes equivalent to 

an effective non-interacting fermionic Hamiltonian. with 
parameters bo, X (where &o = {b^) = (b)) which are deter- 
mined by minimization of the free energy. Thus, we can 
use the two-terminal Landauer formula to calculate the 
conductance, giving at zero temperature 0] 
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where Go is 2e^/ft, and 
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where = = 5ij{uj - E^^ - A) Ms the prop- 

agator of the j-th state of 3 particles with spin pro- 
jection (T in the isolated ring, and the functions A'^^, 

L^iyico - El - A) for 



are A^,, = 1 - b^V^gl^^ 



77 = 77 



and A- = a^^a^^ Jicj - E^ - A) 



for rj ^ 77', where g^\uj) is the Green's function at site 
77 of the corresponding isolated lead = 0). 

In Fig. [2] we show G as a function of magnetic flux in 
the non-magnetic regime E^^"^ < E^^J , for different values 
of V and without electric field (a = 0) . The conductance 
is even with flux 21| and therefore, it is enough to show 
G in the interval < (j) < (j)o/2 (or < <^ab < tt). 
In this regime correlations play a minor role and one 
expects that the JB formula [7|, which is exact in the 
non-interacting case, gives accurate values for G. Our 
results show the same qualitative behavior. In fact, for 
small V it can be demonstrated that both approaches are 



equivalent in this regime. For 
to destructive interference. 



< 



)/2, G vanishes due 



The difference E, 



(4) 



EqJ can be reduced and turned 



negative applying a negative gate voltage. The most 
important results of this work are those obtained in 
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Figure 2: Conductance as a function of magnetic flux for 
a — 0, Vg — 0.8t, U = 2t and several values of V. Full lines: 
our formalism. Dashed lines: JB expression. 



this case i.e., when the ring is in the mixed valence 
or Kondo regime. Results for a = are presented in 
Fig. [3] (a). For small enough A/\E^^^ 
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Figure 3: (a) Conductance and (b) occupancy as a function 
of flux for ci = 0, Vg=0, C/ = 6f and several values of V . 

), charge fluctuations are 
frozen and a clear signature of Kondo physics is displayed 
in the characteristic plateau in G at the ideal conduc- 
tance Go (the unitary limit) 01 ■ This is shown in the 
figure for the smaller values of V at small fluxes. The de- 
pendence of G with flux, is related with the correspond- 
ing dependence of the energy levels and matrix elements 
with (\). For larger V and ^ab ^ ti", the system is in the 
intermediate valence regime, as reflected in Fig. [3] (b) 
in which the total occupancy of the configuration with 
three particles = J2jcrif^f<^) shown. Therefore, the 
conductance deviates from the unitary limit. 

Independently of the other parameters, G vanishes at 
^AB = TT. Within our formulation, at this point the 
states of the 77, = 3 configuration become doubly degen- 
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erate between states of different parity. The matrix ele- 
ments Oj^ entering Eq. ([9|) have the same modulus but 
different sign, therefore producing a complete destructive 
interference inside the absolute value. To our knowledge, 
there are no calculations so far showing at the same time 
this destructive interference and ideal conductance in the 
Kondo regime. The JB expression gives values below 
O.lGo for all ^ab and parameters of Fig. [3l Previous 
mappings to the Anderson model displayed the Kondo 
physics, but did not capture the interference P. [l^. 




Figure 4: Conductance as a function of flux for V — 3.5t and 
difi'erent values of a. Other parameters as in Fig. [3] 

The effect of the SOC on the total conductance is dra- 
matic in the Kondo regime. The results presented in Fig. 
m show dips (additional to that of ^ab = tt) which are 
larger as a grows. The main difference with the case 
a = is that n/| rifi therefore producing a partial de- 
struction of the Kondo resonance, mimicking the effect 
of a Zeeman term. This effect is larger for lower A (when 
the system is deeper inside the Kondo regime), which for 
the parameters of Fig. IH corresponds to ^ab ~ ±0.37r. 
For ^ab = TT, complete cancellation is not achieved. This 
can be interpreted as an effect of quasiparticles acquiring 
relative phases differing slightly from tt due to the effect 
of the additional AC phase (see inset in Fig. [3]). 

Another important effect of the SOC in the Kondo 
regime is that it leads to currents with significant spin 
polarization. If a spin a (up or down) in the quantization 
direction ng is injected in the ring at the right lead ((/? = 
0) it comes out at the left lead {ip = tt) with spin a 
in the direction n'g = (sin0, 0, cos6') or vice versa fl8| . 
The corresponding conductance Go- is spin dependent, 
as shown in Fig. [5l The ratio of the conductances can 
reach a factor 2 or larger with ideal G| (G\ ) for flux 
^AB = O.Wtt (-O.IStt) and rather smaU a [2l[. For 
these values, the z component of the quantization axis 
for any ip is larger than 0.99 

In summary, we have presented an approach to calcu- 
late the conductance through a ring of interacting QDs 
threaded by a magnetic flux and with spin-orbit cou- 
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Figure 5: Conductance for each spin as a function of flux for 
the same parameters as in Fig. |4]and different values of a. 



pling a in the Kondo, mixcd-valcncc and non-magnetic 
regimes. The effects of a are incorporated into Aharonov- 
Casher phases using a gauge transformation that leads 
to the Hubbard Hamiltonian Eqs. (|4]) and ([5]). Using a 
method based on a mapping of the relevant exact eigen- 
states of the ring onto an effective multilevel Anderson 
impurity and with the use of a slave-boson representa- 
tion in the saddle-point approximation, we are able to 
describe the properties of the ring connected to the leads. 
The method is valid for small values of the coupling be- 
tween rings and leads V and small values of magnetic 
field B, such that the Zeeman energy is much less than 
Tk- When the ring is in the Kondo regime, we obtain 
ideal conductance for a = and magnetic flux far from 
half a flux quantum, for which there is complete destruc- 
tive interference. The effect of a small non-vanishing a 
is to produce a progressive destruction of the Kondo ef- 
fect, decreasing the conductance and leading to a strong 
spin dependence of it. Extensions to include the Zeeman 
term or other interacting systems with local interactions 
are straightforward. 
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